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Summary 


A detailed derivation of the equations governing the pressure in a generic liquid-fuel launch 
vehicle tank subjected to uniformly accelerated motion is presented. The equations obtained are 
then for the Space Shuttle Superlightweight Liquid-Oxygen Tank at approximately 70 seconds 
into flight. This generic derivation is applicable to any fuel tank in the form of a surface of 
revolution and should be useful in the design of future launch vehicles. 


Introduction 

An important consideration in the design of launch vehicles is the pressure exerted on the 
wall of a fuel tank during flight. For example, studies of the Space Shuttle Superlightweight 
Liquid-Oxygen (L0 2 ) Tank shown in figure 1 have been presented that include the effects of the 
L0 2 pressure on its nonlinear and buckling behavior under flight loads. 1 2 In these very complex 
studies, the LO, pressure acting on the tank wall was approximated by using several simplifying 
presumptions that are based on the specific characteristics of the vehicle dynamics, at specific 
times in the flight trajectory. For example, stratification of the cryogenic fluid was neglected and 
the fluid was presumed to be homogeneous and, as a result, to have a uniform density. 
Additionally, fluid sloshing effects were neglected. The fluid was taken to be incompressible and 
inviscid, and flow effects were neglected. Moreover, every particle of fluid in the tank was 
presumed to be subjected to uniformly accelerated motion. 

The objective of this memorandum is to present a detailed derivation of the fomulas used in 
references 2 and 3 to simulate the the LO, pressure field in the Space Shuttle Superlightweight 
Liquid-Oxygen Tank. Toward this objective, a generic presentation of the fundamental equations 
is given first, followed by an example problem that corresponds to approximately 70 seconds into 
flight, given in reference 1 . This generic derivation is applicable to any fuel tank in the form of a 
surface of revolution and should be useful in the design of future launch vehicles. 


Derivation 

In the analysis presented subsequently, the fluid is presumed to be incompressible and 
inviscid, and to have a uniform density. The tank and the fluid are presumed to be moving together 
as a uniformly accelerated rigid body without sloshing, and the tank has an internal ullage 
pressure p u . For these presumptions, the pressure field, p(x, y, z), is given by (see reference 3) 

-Vp + pg = pa (1) 

In this equation, p is the density of the fluid, a is the acceleration vector, pg is the body force 

associated with gravity, and V is the gradient operator for the (x, y, z) coordinate system fixed to 
the generic tank shown in figure 2. Associated with the (x, y, z) coordinates are the corresponding 
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unit vectors (i ?£\ • It is convenient to express equation (1) as 

V ’ J ’ / 


- Vp = pa 


(2a) 


where 


a = a - g 


(2b) 


is the net acceleration vector. 

To obtain the desired expression for the pressure field acting on the tank wall, it is convenient 
to first determine the magnitude and direction of net lateral acceleration, denoted by r|, as shown 
in figure 3. The magnitude of the net lateral acceleration is given by 

= \A*y + «z (3) 

and the corresponding direction is given by the unit vector 


fj = cost); j + simp k 


(4a) 


where 


cosrp = 



(4b) 


and 


simp = 


Oz 

a„ 


(4c) 


Next, the coordinate system (x, t], 5) is introduced, as shown in figure 3. The (0, p, £) coordinates 
shown in this figure span the same plane as the (0, y, z) coordinates. With the (x, p, t) coordinates, 
equation (1) is expressed in a simpler component form as 


dp 

dx 


= PW 


(5a) 


dp 

dp 


= P« n 


^ = 0 


(5b) 

(5c) 


From equation (5c), it follows that the functional dependence of the pressure field is given by 
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P = P(X,T1) 

On a surface of constant pressure, the pressure differential dp must obey 

dp = ^E dx + iE d „ + ^Edi; = 0 

ax a,| 1 dt, 

Using equation (5c) with equation (7) yields 

dp 

dx _ _ dq 
dr| dp 

dx 

Substituting equations (5a) and (5b) into equation (8) yields 

dx _ a ri 
dri a x 

which is easily integrated to yield the equation of plane 


( 6 ) 


(7) 


( 8 ) 


(9) 


x + ^q = C (10) 

<*x 

where C is a constant. It is useful to point out that the net acceleration vector in the x-q plane is 
given by 

a = a x i + a, 1 q (11) 

and its slope in this plane is tx x /a , which is the negative reciprocal of the slope of the plane of 
constant pressure given by equation (10). Thus, the planes of constant pressure are perpendicular 
to the net acceleration vector in the x-q plane. The angle of inclination, (3, of the lines of constant 
pressure is given by 

tan P = - ^ (12) 

as seen in figure 4. By using equations (3) and (9), equation (12) becomes 
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In addition, figure 4 indicates that 

cos(3 = 

a 


and 


sin|3 - — 7 - 
a 


(14a) 


(14b) 


To obtain an expression for the pressure gradient that can be easily integrated, it is convenient 
to let (x, rj) denote the coordinates of points in the plane perpendicular to the planes of constant 
pressure, as shown in the figure 4. For this coordinate system, equation (2a) is expressed simply as 



(15) 


Integrating this equation gives 

p = - p| 3 |x + C 

where C is a constant. At the top of the fluid surface, x = x f , where 


(16) 


x f = x f cos|3 — x f 


a, 

131 


(17) 


where equation (14a) has been used. At x = x f , the pressure field is equal to the ullage pressure 
p u , as depicted in figure 4. Therefore, equations (16) and (17) yield 


(18a) 


_► 

( a x 

a 

m _x ) 


for 


X < Xf = Xf 


a, 

lal 


(18b) 


The coordinate x is expressed in terms of the (x, q, £,) coordinates by 


x = cos|3 x + sin|3 q - x + — q 

a a 


Substituting equation (19) into equations (18) gives 


( 19 ) 
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a 

P = Pu + P a x| Xf - X L T] 

a v 


(20a) 


for 


^ a ’i 

x H L ri < x f 

a. 


(20b) 


Next, the r| coordinate is expressed as (see figure 3) 


a, a. 


r| = cosip y + simp z = — y H z 




( 21 ) 


where equations (4) have been used. Substituting equation (21) into equations (20) gives 


a a z 

P - Pu + P a x| x f - x y z 

a. a v 


(22a) 


for 


a y a. 


xh y i z<x f 

a, a. 


(22b) 


For a tank in the form of a surface of revolution, the y and z coordinates of a point (x, y, z) 
on the surface are expressed as 


y = r(x)cos0 

(23a) 

z = r(x)sin0 

(23b) 


where r(x) is the horizontal radius of the tank wall and 0 is the polar angle, as shown in figure 
5. By using equations (23) with equations (22), the pressure field acting on the tank wall is given 
by 

p(x, 0) = p u for x>x s (x, 0) (24a) 

and 


where 


p(x, 0) = p u + pa x x s (x, 0) - x] for x < x s (x, 0) 


(24b) 


x,(x, 0) = Xf - r(x) 


„ Ww 

— cos 0 h sin 0 

a x a x 


(24c) 
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Example Problem 


Formulas similar to equations (24) were used in the structural analysis of the Space Shuttle 
Superlightweight External Tank (SLWT). In particular, an early booster ascent loading condition 
was examined in reference 1 that corresponds to approximately 70 seconds into flight and a LO, 
tank that is approximately seven-eighths full. At this instant of the trajectory, the acceleration field 
was approximated by a uniform field that is given in terms of the coordinates shown in figure 5 
by a x = 2.01 lg, a = 0.049g, and a z = 0.440g, where g is the magnitude of uniform gravitational 
acceleration. For this acceleration field, equations (3) and (4) give rf = 83.65 degrees for the plane 
of the resultant lateral acceleration vector. Likewise, equation (13) gives (3 = 12.42 degrees for 
the inclination of the fluid surface, in the x-r| plane, shown in figure 4. These two angles are 
independent of the shape of the tank. 

The SLWT LO, tank, without its foam insulation, is shown in figure 6 and the corresponding 
geometry and dimensions of a corresponding idealized wall reference surface are shown in figure 
7. The major tank-wall components consists of a forward ogive, an aft ogive, a cylindrical barrel 
section, and an aft dome. The aft dome consists of a truncated elliptical dome section with a 
spherical cap. Each of these tank elements has a circular cross section. 

Determination of the pressure acting on tank wall requires an expression for the horizontal 
radius r(x), appearing in equation (24b), for each tank component. The geometry of the ogive 
section, consisting of the forward and aft ogives, is shown in figure 8. The ogive is formed by 
rotating a segment of an eccentric circle about the x-axis. This tangent line of the segment at 
x = h 2 is parallel to the x-axis. Inspection of this figure reveals 

x = h 2 + Rsincj) (25a) 

and 

r(x) = R(coscl) - 1) + r 2 (25b) 


From equation (25a) it follows that 


COS(|) 


1 

R 


V R - ( x - h 2 ) 


(26) 


Substituting equation (26) into equation (25b) gives the desired expression 

r(x) = a^R - (x - h 2 ) + r 2 - R for h 2 <x<h,+h 2 
For a given value of r, shown in figure 7, the height h, is found from equation (27) as 

h, = <J(r 2 - r,)[2R - r 2 + r,] 


(27) 


(28) 
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For the barrel section, inspection of figure 7 gives 


r(x) = r, for 0<x<h, (29) 

The corresponding expression for the truncated elliptical section of the aft dome is obtained 
directly by using the equation for an ellipse in x-r coordinates along with the geometry given in 
figure 7. The result is 



for - h 3 < x < 0 


(30) 


Likewise, the expression for the spherical part of the aft dome is obtained by using the equation 
for a circle whose origin is at r = 0 and x = R - h 3 - h 4 , where R is the radius of the spherical 
cap shown in figure 7. From this equation, it follows that 


r(x) = y R( - (x - R s + h 3 + h 4 ) for -(h 3 + h 4 ) 


< x < 


(31a) 


where 


R s = 


2 i 2 

£3 + h 4 
2h 4 


(31b) 


A "to scale" graphical image based on equations (27)-(31) is shown in figure 9. For this image, a 
value of Tj = 0 was used to obtain the sharp point at the nose of the tank. The corresponding value 
of h, given by equation (28) is 418.5 in. 

The approximate pressure distribution on the interior surface of the LO , tank is given by 
equations (24), where x is the axial coordinate shown in figure 7 and r(x) is given by equations 
(27)-(3 1). The fill level for this example is given by x = x f = 305.9 in., and the density p is given 
by pg = y, where y = 0.04123 lb/in 3 is the specific weight of the LO,. The ullage pressure for this 
example is p u = 19.6 psi and the corresponding tank dimensions are given in figure 7. Pressure 
calculations were made for these specific values and for a closed, pointed ogive with r, = 0, and 
are shown in figures 10 and 11. The results in figure 10 show the variation of the pressure acting 
on the shell wall as a function of the polar angle 0. Five curves are shown in the figure that 
correspond to different axial positions along the tank. The straight line shown for x = 500 in. is 
in the unfilled region of the ogive where the pressure is equal to the 19.6 psi ullage pressure, and 
is independent of 0. The curves shown for x = 250, 100, and 0 in. correspond to a location in the 
ogive below the fluid surface, a location within the barrel section, and the location where the 
barrel and aft dome meet, respectively. The sinusoidal shape of these curves corresponds to the 
inclination of the fluid surface caused by the uniformly accelerated motion. The straight line 
shown for x = -1 14 in. is almost at bottom of the aft dome, where the variations of the pressure 
with 0 are negligible. 
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The results in figure 1 1 show the variation of the pressure acting on the shell wall as a 
function of the axial coordinate x. Two curves are shown in this figure that correspond to 
0 = 83.65 and 263.65 degrees. These two tank meridians span the plane of maximum lateral 
acceleration. The differences in the two curves are also associated with the inclination of the fluid 
surface. The two curve intersect at x = -1 14.9 in., the center of the aft dome, where a maximum 
pressure equal to 54.5 psi occurs. 


Concluding Remarks 

A detailed derivation of the equations governing the pressure in a generic liquid-fuel launch 
vehicle tank subjected to uniformly accelerated motion has been presented. The equations 
obtained have been applied to detemine the pressure acting on the shell wall of the Space Shuttle 
Superlightweight Liquid-Oxygen Tank at approximately 70 seconds into flight. This generic 
derivation is applicable to any fuel tank in the form of a surface of revolution and should be useful 
in the design of future launch vehicles. 
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Figure 1. Space Shuttle external tank components. 
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Figure 2. Baseline coordinate frame of generic fuel tank. 



Figure 3. Coordinate frame of generic fuel tank associated with the resultant 
lateral acceleration. 
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Figure 4. Coordinate frame of generic fuel tank associated with lines of constant 
pressure in the plane of resultant lateral acceleration. 



Figure 5. Cylindrical coordinate frame for a generic fuel tank in the form of a 
surface of revolution. 
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Figure 6. Space Shuttle Superlightweight Liquid-Oxygen Tank. 
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Figure 7. Geometry and dimensions of Space Shuttle Superlightweight Liquid-Oxygen Tank. 
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X 



Figure 8. Geometry of the ogive. 
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Figure 9. Graphical output of equations (27)-(31) for the Space Shuttle Superlightweight 
Liquid-Oxygen Tank with the dimensions given in figure 7, except r, = 0. 



Figure 10. Pressure acting on the shell wall of the Space Shuttle Superlightweight 
Liquid-Oxygen Tank as a function of the polar angle 0 (see figure 5). 
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Axial coordinate, x, in. 


Figure 11. Pressure acting on the shell wall of the Space Shuttle Superlightweight 
Liquid-Oxygen Tank as a function of the axial coordinate x (see figure 5). 
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